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Abstract. For applications in computing, Bezier curves are perva- 
sive and are defined by a piecewise linear curve £ which is embedded 
in R"^ and yields a smooth polynomial curve C embedded in R"^. It is 
of interest to understand when C and C have the same embeddings. 
One class of counterexamples is shown for C being unknotted, while 
C is knotted. Another class of counterexamples is created where C is 
equilateral and simple, while C is self-intersecting. These counterex- 
amples were discovered using curve visualizing software and numerical 
algorithms that produce general procedures to create more examples. 



Computational topology lies properly within the broad scope of applied gen- 
eral topology and depends upon a novel integration of the 'pure mathematics' 
of topology with the 'applied mathematics' of numerical analysis. Compu- 
tational topology also blends general topology, geometric topology and knot 
theory with computer visualization and graphics, as presented here. For the 
computational cases considered here, ideas from general topology, geometric 
topology and knot theory are complemented by numeric arguments in a novel 
integration of the 'pure' field of topology with the 'applied' focus of numeri- 
cal analysis. Additionally, aspects of computer visualization and graphics are 
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(a) Unknotted C with knotted C 



(b) Zoomed-in view of C 



Figure 1. Unknotted C with knotted C 



1. Introduction 



(a) Equilateral, simple C with self-intersecting C (b) Zoomed-in self-intersection 



Figure 2. Equilateral, simple C with self-intersecting C 

incorporated into the proofs. Much of this work was motivated by modeling 
biological molecules, such as proteins, as knots for visualization synchronized 
to simulations of the writhing of these molecules. 

Attention is restricted to when C is closed, implying that £ is also closed. 
As C is created from £, it is natural to ask which topological characteristics 
are shared by these two curves, particularly as the control polygon often serves 
as an approximation to the Bezier curve in many practical applications. How- 
ever, topological differences between a Bezier curve and its control polygon 
can exist and it is natural to develop counterexamples to show these topologi- 
cal differences. These counterexamples extend beyond related results, while we 
introduce new computational methods to generate additional counterexamples. 

In Section |2j we present a counterexample of Bezier curve and its control 
polygon being homeomorphic, yet not ambient isotopic. To develop this coun- 
terexample, we created and used a computer visualization tool to study topo- 
logical relationships between a Bezier curve and its control polygon. We viewed 
the images to motivate formal proofs, which also rely on numerical analysis and 
geometric arguments. 

In Section [3) we present numerical techniques to create a class of topological 
counterexamples - where a Bezier curve and its control polygon are not even 
homeomorphic, as the Bezier curve is self-intersecting while the control polygon 
is simple. We exhibit self-intersection by a numerical method, which finds 
the roots of a system of equation. We freely admit that these roots are not 
determined with infinite precision, but such calculations on polynomials of 
degree 6, as in these examples, typically elude precise calculation. We argue 
that two primary values of our method for these approximated solutions are 

(1) as a catalyst to alternative examples that may admit infinite precision 
calculation and rigorous topological proofs, and 

(2) as having specified digits of accuracy - typically crucial for acceptable 
approximations in computational mathematics and computer graphics. 

Using the visualization tool described in Section [4j we viewed many examples 
where the Bezier curve was simple, while its control polygon was equilateral 
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and simple. It is well known that Bezier curve can be self-intersecting even 
when its control polygon is simple, but we conjectured that the added equilat- 
eral hypothesis would imply that both curves were simple. While this visual 
evidence was suggestive, we present a general numerical approach in Section |3] 
that supports a contrary interpretation. This example provides guidance for 
designing appropriate approximation algorithms for computer graphics. 

1.1. Mathematical Definitions. The definitions presented are restricted to 
R^, as sufficient for the purposes of this paper, but the interested reader can 
find appropriate generalizations in published literature. 

Definition 1.1. A knot [4 is a curve in which is homeomorphic to a circle. 

Knots are often described by a knot diagram [12 , which is a planar projection 
of a knot. Self-intersections in the knot diagram correspond to crossings in 
the knot, where each crossing has one arc that is an undercrossing and an 
overcrossing^ relative to the direction of projection. 

Homeomorphism is an equivalence relation over point sets and does not 
distinguish between different embeddings, while ambient isotopy is a stronger 
equivalence relation which is fundamental for classification of knots. 

Definition 1.2. Let X and Y be two subspaces o/R^. A continuous function 

H X [0, 1] ^ R^ 
is an ambient isotopy between X and Y if H satisfies the following: 

(1) i^(-,0) is the identity, 

(2) hIx, 1) = Y, and 

(3) \/t G [0, l],i^(-,t) is a homeomorphism from R^ onto R^. 
The sets X and Y are then said to be ambient isotopic. 

Definition 1.3. Denote C{t) as the parameterized Bezier curve of degree n 
with control points Pm G R^^ defined by 

n 

C{t) = Y,B^At)P^^t^ [0,1] 

i=0 

where 5i,n(t) = t'{l - tY''. 

1.2. Related Literature. A Bezier curve and its control polygon may have 
substantial topological differences. It is well known that a Bezier curve and its 
control polygon are not necessarily homeomorphic [14 . Recently, it was shown 
that there exists an unknotted Bezier curve with a knotted control polygon 
|6]. A specific dual example has also been shown [16 of a knotted Bezier 
curve with an unknotted control polygon. However, the methodology was a 
visual construction without formal proof and is not easily generalized. Software, 
Knot-Splinc-Vis^ developed by authors Marsh and Peters, was used to find 
another example, where the methodology can more easily be generalized and 
Knot-Splinc- Vis is publicly available [13J . 
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Much of the motivation for considering these counterexamples came from 
apphcations in scientific visuahzation [9l [TOj [11] . A primary focus was to es- 
tabhsh tubular neighborhoods of knotted curves so that piecewise linear (PL) 
approximations of those curves within those neighborhoods remained ambient 
isotopic to the original curves. This was initially considered for approximations 
used in producing static images of these curves, but it became readily apparent 
that these same neighborhoods also provided bounds within which many per- 
turbations of these models remained ambient isotopic under these movements. 
That theory is being applied to dynamic visualizations of molecular simula- 
tions, where the neighborhood boundaries permit convenient warnings as to an 
impending self-intersection, as of possible interest to biologists and chemists 
who are running the simulations. 

Previous work [15 in knot visualization provides a rich set of data for PL 
knots, where each edge is of the same length. The interface of KnotSpUne-Vis 
was designed to import this equilateral PL knot data and then generate the 
associated Bezier curves. This matured into an empirical study of dozens of 
cases, where all examples examined yielded simple Bezier curves for these sim- 
ple equilateral control polygons. This raised the question of whether the pres- 
ence of equilateral edges in the control polygon might be a sufficient additional 
hypothesis to ensure homeomorphic equivalence with the Bezier curve, as the 
previously cited examples [14 did not have equilateral control polygons. 

2. Unknotted C with knotted C 

In order to produce a knotted Bezier curve with an unknotted control poly- 
gon, we invoked KnotSpline-Vis with an example (Figure [3| of an unknotted 
Bezier curve, where the total curvature appeared to be larger than An. (The to- 
tal curvature being larger than An is a necessary condition of knottedness.) We 




(a) Unknotted C & The V (b) Zoomed-in View of C 



Figure 3. Visual experiments 
experimented on this example by moving control points to construct a Bezier 
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curve that visually appeared to be knotted. We then moved control points to 
unknot the control polygon while keeping the Bezier curve knotted. In th e end 



we obtained a Bezier curve and the control polygon (of degree 10) (Figure 1(a) ) 
defined by the control points {Pq, Pi, • • • , Pg, Pq} listed below: 

(-5.9, 4.7, -6.2), (10.3, -1.1, 8.9), (-2.6, -12.4, -6.3), (-10, 7, -0.3), 
(1.9, -12, -0.6), (11.2, 7.5, -7.6), (-15.3, -1.7, -4.1), (-11.7, 20, 3.5), 
(17.9, -1.1, 2.9), (2.9, -13.7, 4.8), (-5.9, 4.7, -6.2). 

The 3D visualization offers only suggestive evidence that the above Bezier 
curve is knotted while the control polygon is unknotted. We provide rigorous 



mathematical proofs of these properties in Sections 2A_ and |2.2| Generally, 
proving knottedness or unknottedness can be difficult [8 , but is accessible for 
the counterexample here. 

2.1. Proofs of the Bezier curve being knotted. We prove that the Bezier 
curve is a trefoiQ We orthogonally project the curve onto x-y plane. We 
then shown that there are three self-intersections in the projection and these 
intersections are alternating crossings in 3D. Since projections preserve self- 
intersections, this curve can have no more than 3 self-intersections, but these 
self-intersections in x-y plane are shown to have pre-images that are 3D cross- 
ings, so the original curve has no self- intersections. 




Figure 4. The 2D projection of the knot 

Since the 2D curve in Figure [4] has degree 10, we use a numerical method im- 
plemented by MATLAB function 'fminsearch' to find the parameters where the 



■"^A trefoil is a knot with three alternating crossings (12). 
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curve intersects with itself. We provide the numerical codes in Appendix ] A. 1[ 
and we provide the data used to find these parameters in Appendix |A.3[ The 
pairs of parameters (labeled in order as ti, . . . ,^6) of the self-intersections are 
listed below: 



[h = 0.0306,^4 = 0.5573], [^2 = 0.1573,^5 = 0.9244], [ts = 0.3731,^6 = 0.9493]. 

Next we prove that these 2D intersections are projections from three alter- 
nating crossings in 3 D. The above parameters are substituted into the Bezier 
curve (Definitio n|1.3|) to get pairs of points (numerical codes for this calculation 
are in Appendix |A.2p : 

[C{ti) = (-2.0539, 2.8001, -2.6929), C(t4) = (-2.0530,2.7987,-2.0143)], 

[C{t2) = (0.4376, -2.5212, -0.0576), ^(ts) = (-0.4364, -2.5206, -0.5547)], 

[C(t3) = (-1.3613, -1.4239, -2.2944), ^(te) = (-1.3624, -1.4232, -1.9067)]. 

The alternating crossings follow from comparing the z-coordinates in each 
pair. Precisely, according to the parameters given above, the tracing of the six 
points in order is 

C(ti),C(t2),C(t3),C(t4),C(t5),C(t6), 

and the crossings at these points are 

under ^ over, under, over, under, over. 

2.2. Proof of the control polygon being unknotted. To prove that the 
control polygon P = (Pq, Pi, • • • , Pg, Pq) is an unknot, it is necessary to show 
that P is simple and unknotted. We directly tested each pair of segments 
of P for non-self-intersection and those calculations can be repeated by the 
interested reader. We prove unknottedness using a 3D push, as the obvious 
generalization of a 2D function from low- dimensional geometric topology [5]. 
We restrict attention to a median push, as defined below. The full sequence of 
5 median pushes is explicated, where the first 4 median pushes are equivalently 
described by Reidemeister moves [12 . 

Definition 2.1. Suppose APj_iPjPj+i is a triangle determined by non-collinear 
vertices Pj-i,Pj and Pj+i of P. Push the vertex Pj, along the corresponding 
median of the triangle to the middle point of the side Pj_iPj+i. We call this 
function a median push. 

We depict the sequence of pushes used in Figure [5j showing the 3D graphs 
of the PL space curves after the median pushes. These graphs show, at each 
step, which vertex is pushed and its image. For example, in Figure |5(a)[ the 
vertex P3 is pushed to P3 and the resultant polygon after this push is shown by 
Figure [5(b)| Figures |5 (a) ||5 (b) [5 (c) | and |5 (d) | have corresponding Reidemeister 
moves, as those pushes eliminate at least one crossing. Using the published 



notation [12 , Figure 5(a) depicts a Reidemeister move of Type 2b. Similarly, 
Figure [5 (b)| depicts a Reidemeister move of Type lb; Figure [5(c) has Type 2b 
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(a) The initial control polygon 
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(c) After the 2nd push 
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(b) After the 1st push 




(d) After the 3rd push 




(e) After the 4th push (f) The unknot after these pushes 

Figure 5. Pushes and PL curves in 3D 



and Figure 5(d) has a move of Type 2b, fohowed by a move of Type lb. The 
final push to achieve Figure [5(f)] does not correspond to any Reidemeister move 
as no crossings are changed, but it is included to have a polygon with only five 
edges, which necessarily must be the unknot [1 . 

We prove that P is unknotted by showing that P is ambient isotopic to 
the unknotted PL curve shown in Figure 5(f) As a sufficient condition, we 
show that the pushes do not cause intersections [?. (We gained significant 
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intuition for specifying the sequence of pushes by visual verification with our 
3D graphics software capabihties to translate, rotate and zoom the images.) 
We now present th e form al arguments. 



Consider Figure 5(a) , where P3 is pushed to P3 (the middle poinlnof P2P4)' 



We show that any segments other than P2P3 and P3P4 of P do not intersect 
the triangle AP2P3P4 or the triangle interior. 
We parameterize each segment by: 



FJR^i : P, + (P,+i - Pi)t, t e [0, 1] 
for z = 0, 1, • • • ,9 and let Pio = Pq. Then the points given by 

P3 + a{P2-P3) + b{P4-P3), 

for a, 6 > and a -\- b < 1 are on the AP2P3P4 and contained in its interior. 
Hence Pi + (Pi+i — Pi)t intersects AP2P3P4 or its interior if and only if 

(2.1) P, + (P,+i - P,)t = P3 + a(P2 - P3) + 6(P4 - P3) 

has a solution for some t G [0, 1] and a, 6 > and a -\- b < 1. 



For each i = 0,1,-- - ,9, we solve Equation 2.1 (a system of 3 linear equa- 
tions) for a, b and t with the above constraints. Calculations show there is 
no solution for each system, and hence Pi + (Pi+i — Pi)t does not intersect 
AP2P3P4 or its interior for each z = 0,l,---,9. Thus it follows that the push 
does not cause any intersections. 

Similar computations verify that the other pushes do not cause any inter- 
sections, thus establishing the ambient isotopy. 

3. Equilateral, simple C with self-intersecting C 

We visualized many cases of simple, closed equilateral control polygon^ that 
all appeared to have unknotted Bezier curves. Many of these control polygons 
were nontrivial knots. Prompted by this visual evidence, we conjectured that 
any simple, closed equilateral control produces a simple Bezier curve, where 
some examples are shown in Figure |6] 

We now present numerical evidence to the contrary. As noted in Section [TJ 
the degree 6 polynomials make precise computation difficult, so we do not 
provide a completely formal proof, independent of numerical methods. A le- 
gitimate concern is whether the numerical approximation produces coordinates 
for a 'near' intersection, subject to the accuracy of the ffoating point compu- 
tation. We cannot refute that possibility, but we argue that in the context of 
graphical images, the level of approximation produced is often sufficient. In 
particular, we have parameters in the code to adjust the number of digits of 
accuracy. This is level of user-deffned precision is often accepted as sufficient 
for visualization |9j. The user can then set the graphical resolution so that 
points that are determined to be the same within some acceptable numerical 
tolerance will also appear within the same pixel. 



9 ^ 

It is not necessary to push P3 to the middle point. Any point along P2P4 would suffice. 
"^Throughout this section, all control polygons are simple. 



Computational Topology Counterexamples 



9 



Figure 6. Unknotted Bezier curves with equilateral control polygons 

3.1. Intuitive overview. We create many examples to test and retain only 
those that satisfy all three criteria listed in Section |3.4[ We begin the cre- 
ation of a closed equilateral polygon by setting Pq = (O^O^O)- We then take 



{Qo^Qi^ ' ' ' 7<7n-i} (Equation 3.3) from the unit sphere so that 



Pl=Po+gO, P2=Pl^qU Pn = Pn-l ^ qn-l^ 



We ensure that the polygon is closed in Section 3.4 item (3). 



We consider a sufficient and necessary condition for a Bezier curve being self- 



intersecting (Equation 3.2). Since we want the equilateral polygon to define a 
Bezier curve that is self-intersecting, we not only select {qo^ Qir ' ' •> <?n-i} from 
the unit sphere as above, but also select them such that Equation |3.2| is zero 
for some parameters s and t. Consequently the set of control points generated 
determines a closed equilateral control polygon and a self-intersecting Bezier 
curve. 

3.2. Necessary and sufficient condition for self-intersection. We rely 
upon the following published equation [3 for necessary and sufficient conditions 
for self-intersection of a Bezier curve 

with the domain D = {{s,t) : s < > 0,{s,t) ^ (0,0)}. A Bezier 

curve defined by C{t) is self-intersecting if and only if there exist s and t in the 
domain D such that S{s^t) = 0, , with an alternative formulatiorj^ given by 



^One needs to write out the 3 Equation (6)] to obtain Equation 3.2 
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(3.2) 

^ i=o j=o V ^ / fe=o 

where 

(3.3) qi=P,^i-P, /or z = {0,1,- ,n-l}. 

3.3. A representative counterexample generated. We present a single 
numerical counterexample, noting that while only one counterexample is pre- 
sented, the numerical algorithm implemented can be used to find many such 
examples. We list six distinct control points (the seventh control point is equal 
to the first control point) that determine an equilateral simple C and a self- 
intersecting C (shown in Firgure|2|, as generated by the algorithm described in 
detail in Section l3^ 



(0, 0, 0), (0.0305, 0.0810, 0.9962), (-0.2074, -0.2671, 1.9030), 
(-0.1792, -0.3402, 0.9063), (0.0189, 0.0782, 0.0185), (0.1557, 0.2329, -0.9600). 

We verify that C is simple by considering all pairwise intersections of the 
segments of this control polygon. The self-intersection of the Bezier curve 
occurs at a point that is numerically approximated as 

[s,t] = [0.2969,0.0633] 

where correspondingly, 

S{s,t) = (-0.0003861,-0.000097,0.0001462) ^ (0,0,0). 
The error occurs because of numerical round off on s, t and the control points. 

3.4. The numerical method for generating counterexamples. We pro- 
vide the numerical codes and data used in Appendix [B] Given a control poly- 
gon, we can determine whether a self-intersection of the Bezier curve occurs 



by determining whether S (Equation 3.2) has a root in D. We consider as a 



function S{s^t^q) where q = {qi}^^Q, so that finding a self-intersecting Bezier 
curve with an equilateral control polygon is equivalent to determining s^t and 
q such that the following are satisfied: 

(1) S{s,t,q) = where s,t e D; 

(2) \\qi\\ = r for each z G {0, 1, • • • , n — 1}, where || • || is the Euclidean 
norm; 

(3) "^27=0 ~ ^7=0 — Pi = ^ since P needs to be closed. 

We assume r = 1 without loss of generality since the value of r can be 
adjusted by scaling. Throughout the provided codes, n is always the degree of 



the Bezier curve. We give the code for function S in Appendix B.l , where [5, t] 
is labeled as u and q as [a, 6, c]. 



Computational Topology Counterexamples 



11 



The function SF (Appendix B.2 ) takes parameters for s, t and q and outputs 
a floating point value. It is designed to be zero if and only if the above three 
conditions are satisfied simultaneously. 

Precisely since q should be taken from the unit sphere, SF assigns a, 6, c 
values given by 

a = sin{(j))cos{6); 
b = sin{(j))sin{0)] 
c — cos{(j))^ 

where (j) and represent input parameters. In order to satisfy Condition (3) 

above, qn-i is set equal to — Xir=o^^*- -^^^ ^^^^ ^^Yi lkn-i|| fail to 
be equal to 1. So we include the function F (Appendix |B.2[ ) to determine 
whether = 1- The function is designed to be F = llg'n-iH — 1 such that 

\\qn-i\ \ = 1 if and only if F = 0. 
Symbolically, 

(3.4) 5F=||5|| + ||F||, 



where S is given by Equation 3.2 and F = \ \qn-i\ \ — 1. Having the above three 
conditions satisfied simultaneously is equivalent to finding input values such 
that SF = 0. 

Since SF > 0, the minimum of SF is 0. The function SFminimizer (Ap- 



pendix B.3) uses fminsearch (a numerical method integrated in MATLAB) 
to search for the minimum of SF, while returning this minimum and the cor- 
responding values of s,t and q. The user supplied initial guesses for s,t and 
q greatly infiuence the results, so we assign SFminimizer randomized initial 
values M times so that we get M different minimums for different initial val- 
ues. But no matter which initial values we use, as long as we can get "a" 
minimum of 0, then we get the equilateral control polygon which determines a 
self- intersecting Bezier curve. 

The data of finding the counterexample of Section [33| is included in Appen- 
dix [El 

4. Visualizing Bezier curves & their control polygons 

To study the knot types of a control polygon and its Bezier curve, a knot visu- 
alization tool was developed, called Knot Spline-Vis. The tool KnotSpline-Vis 
takes a PL control polygon as input and displays both the PL curve and the 
associated Bezier curve. For these studies, the input was always restricted to 
be a PL curve of known knot type. The functions of Knot_Spline_Vis were 
designed to permit interactive studies of the topological relationships between 
a Bezier curve and its control polygon. The intent is to use these examples to 
stimulate mathematical conjectures as a prelude to formal proofs. 

Some of the standard graphical manipulation capabilities provided are il- 
lustrated in the following Figure [7| and |8) where the rotation capabilities have 
been particularly useful to develop visual insights into the occurrence of self- 
intersections and crossings as fundamental for the study of knots. An editing 
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Figure 8. Scaling 

The control polygon is an initial PL approximation to its associated Bezier 
curve. Subdivision algorithms [14 produce additional control points to further 
refine the original control polygon, converging in distance to the Bezier curve. 



Figure 10 illustrates performance of a standard subdivision technique, the de 



Casteljau algorithm [7 . Users can specify a subdivision parameter. Figure 10 
shows the case with parameter ^. 

5. Conclusions and Future Work 

We use numerical algorithms and graphical software to generate counterex- 
amples regarding the embeddings in for a Bezier curve and its control poly- 
gon. 
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(a) Initial display (b) Some vertices moved 



Figure 9. Moving control points 




Figure 10. Subdivision 



First, we present a class of counterexamples showing an unknotted control 
polygon for a knotted Bezier curve and provide complete formal proofs of this 
condition. We used a spline visualization tool, KnoLSpline.Vis for easy gen- 
eration of many examples to gain intuitive understanding to formalize these 
proofs. Secondly, we used KnotSpline-Vis in formulating the conjecture that 
any simple, closed equilateral control had a simple Bezier curve. We provide 
contrary numerical evidence regarding this conjecture, as a useful guide for 
many applications in computer graphics and scientific visualization. 

This work also shows the importance of continuing investigation 

• theoretically, into an infinite precision proof for a counterexample of the 
simple, closed equilateral polygon conjecture, where these demands will 
likely require further, substantive mathematical innovation and 

• practically, into a user interface for KnotSpline- Vis to permit interac- 
tive 3D editing, as the current text driven interface was cumbersome. 
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6. Web Posting of Supplemental Materials 
Appendices listed below are posted on this webpage: 



http : / / www . math . uconn . edu/~ j ili/Supplemental-materials . pdf 



Appendix A: Numerical codes for knottedness of C (Figure 1(a)): 



• A.l Codes for searching self-intersections in Figure [4j 

• A. 2 Codes for determining under or over crossings; 

• A. 3 Data for searching intersections in Figure [4j 



Appendix B: Numerical codes for searching the Example of Figure [2(a)| 



• B.l Codes for Equation |3. 2 

• B.2 Codes for Equation 



B.3 Codes for searching the roots of the system of Equations|3.2|and|3.4 



B.4 Data for finding the example shown in Figure 2(a) 
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A. Code: knottedness of Bezier curve of Figure [1(a) 

A.l. Code for self-intersections of Figure |4| %The function C2d{t) de- 
fines the projection of the Bezier curve, 
function [value] = C2d(t) 
n=10; 

P=zeros(2,ll); 

P(l,:)= [-5.9, 10.3, -2.6, -10, 1.9, 11.2, -15.3, -11.7, 17.9, 2.9, -5.9]; 
P(2,:)= [4.7, -1.1, -12.4, 7, -12, 7.5, -1.7, 20, -1.1, -13.7, 4.7]; 

sum=0; 
for i=0:n 

sum = sum + nchoosek{n^ z) * * (1 — t)^n — i) * P(l, i + 1); 
end 

vi = sum; 

sum=0; 
for i=0:n 

sum = sum + nchoosek{n^ z) * * (1 — t)^n — i) * P(2, i + 1); 
end 

V2 = sum; 



value = ['^1, '^2]; 

% Use the function 'fminsearch' to find the minimums of fnS{x)^ the zero 
minimums and parameters where the zeros occur are what we look for. 



function [value] = fnS(x) 

value = norm{C2d{x{l)) - C2(i(x(2)), 2); 

function [value] = Smin(sO,tO) 

Max=optimset('MaxFunEvals',le+19); 

comb=@(x)fnS(x); 

u=[sO,tO]; 

[xval^ fval] = f min search {comb ^u^ Max) 

A. 2. Code for determining crossings. %Below is the function of the Bezier 
curve in 3D (Definition |1.3| ). 
function [value] = C(t) 
n=10; 

P=zeros(3,ll); 

P(l,:)= [-5.9, 10.3, -2.6, -10, 1.9, 11.2, -15.3, -11.7, 17.9, 2.9, -5.9]; 
P(2,:)= [4.7, -1.1, -12.4, 7, -12, 7.5, -1.7, 20, -1.1, -13.7, 4.7]; 
P(3,:)= [-6.2, 8.9, -6.3, -0.3, -0.6, -7.6, -4.1, 3.5, 2.9, 4.8, -6.2]; 



sum=0; 
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for i=0:n 

sum = sum + nchoosek{n^ z) * * (1 — t)^n — i) * P(l, i + 1); 
end 

vi = sum; 

sum=0; 
for i=0:n 

sum = sum + nchoosek{n^ * * (1 — t)^n — i) * P(2, z + 1); 
end 

V2 = sum; 

sum=0; 
for i=0:n 

52xm = sum + nchoosek{n^ z) * * (1 — t)^n — i) * P(3, z + 1); 
end 



value = ['yi,'y2,'^3]; 

A. 3. Data for searching intersections in Figure |4| % The Matlab com- 
mands and corresponding results: 

% The initial input (0.03, 0.55) is figured out by the observation and calcula- 
tions on self-intersections in Figure [4j Similarly for the others below. 
Smin(0.03,0.55) 
xval = 0.0306 0.5573 
fval = 2.9567e-04 



Smin(0.15,0.92) 
xval = 0.1573 0.9244 
fval = 1.5848e-04 



Smin(0.37,0.95) 
xval = 0.3731 0.9493 
fval = 1.4637e-04 



B. Code: generating counterexamples like Figure [2(a)| 

B.l. Code for Equation 3.2 function [value] = S(u,a,b,c) 
sum=0; subsuml=0; subsum2=0; 
for i=0:n-l 
for j=0:n-l-i 

% The use of adding 1 to the vectors a, b and c is required by MATLAB. 
for k=0:i subsum2 = subsum2 + nchoosek{i^ /c) * (1 — ii(2))*^*~^^ * u{2)^ * a{j - 
k^l); 
end 

subsuml = subsuml H- subsum2 * nchoosek{n — 1 — i, j) * u{l)^^~^~'^~^'^ * (1 
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u{l)y; 
end 

sum = sum + subsuml; 
end 

vi — (1/n) * sum; 



sum=0; subsuml=0; subsum2=0; 
for i=0:n-l 
for j=0:n-l-i 
for k=0: i 

subsum2 = subsum2 + nchoosek{i^ /c) * (1 — 'u(2))'^*~^^ * u{2)^ * b{j + /c + 1); 
end 

subsuml = subsuml + subsum2 * nchoosek{n — 1 — j) * u{l)^'^~^~'^~^'^ * (1 — 

^1)^; 

end 

sum — sum + subsuml; 
end 

'^2 = (1/n) * sum; 

sum=0; subsuml=0; subsum2=0; 
for i=0:n-l 
for j=0:n-l-i 

for k=0:i subsum2 = subsum2 + nchoosek{i^ /c) * (1 — ii(2))^*~^^ * u{2)^ * c{j + 

fc+l); 

end 

subsuml = subsuml + subsum2 * nchoosek{n — 1 — j) * * (1 — 

end 

5ixm = 5iim + subsuml; 
end 

= (1/n) * Slim; 
value = [vi,V2, vs]; 



B.2. Code for Equation 3.4 function [value] = SF(x,n) 

q=zeros(3,n); p=zeros(3,n+l); 

a=zeros(l,n);b=zeros(l,n);c=zeros(l,n); 

alpha = zeros{l^ 2 * n — 2); 

forz = l:2*n — 2 

alpha(i)=x(i); 

end 

for i=l:n-l 

a{i) = sin{alpha{i)) * cos{alpha{n — 1 + i)); 
b{i) = sin{alpha{i)) * sin{alpha{n — 1 + i)); 
c(z) = cos{alpha{i)); 
q(:,i) = [a(i),b(i),c(i)]; 
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end 

a(n)=0; b(n)=0; c(n)=0; 
for i=l:n-l 

a(n)=a(n)-a(i); b(n)=b(n)-b(i); c(n)=c(n)-c(i); 

q(:,n) = [a(n),b(n),c(n)]; 

end 

u=zeros(l,2); u = [x{2 * n — 1), x(2 * n)]; 
value=abs(F(alpha,n)) +norm(S(u,a,b,c) ,2) ; 

for i=l:n 
for j=l:i; 

p(:,i+l)=p(:,i+l)+q(:,j); 

end 

end 

P 



function [value] =F (alpha, n) 

for i=l:n-l 

a{i) = sin{alpha{i))^cos{alpha{n — l-\-i));b{i) — sin{alpha{i)) ^ sin{alpha{n — 

1 + z)); c{i) = cos{alpha{i))] 

end 

a(n)=0; b(n)=0; c(n)=0; 
for i=l:n-l 

a(n)=a(n)-a(i); b(n)=b(n)-b(i); c(n)=c(n)-c(i); 
end 

value = abs{a{n)'^ + 6(n)^ + c(n)^ — 1); 



B.3. Codes for solving the system of Equations |3.2| and |3.4[ function 

[value] = SFminimizer(n,M) 

Max=optimset(MaxFunEvals,4e+9); 

comb=@(x)SF(x,n); 

xval=zeros(2*n,M) ; 

fval=zeros(l,M); 

for i=l:M 

phi=unifrnd(0,pi,l,n-l); theta=unifrnd(0,2*pi,l,n-l); 
sO=unifrnd(0,l); tO=unifrnd(0,l-sO); 
xO= [phi,theta,sO,tO] ; 

[xval{:, i), fval{i)] = fminsearch{comb, xO, Max); 

%make sure s -\-t = xval{2 * n — 1, i) + xval{2 * n, i) < 1 so that s^t G D. 

if xval{2 * n — 1, i) + xval{2 * n, i) > 1 or xval{2 * n — 1, + xval{2 * n, z) = 1 
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fval(i)=l; 
else 
end 
end 

%Display the minimums and corresponding 5, t and q. 
xval fval 

%Returns the optimal one. 
[C, /] = min{fval)] 
value=[C,I]; 

B.4. Data for finding the example shown in Figure [2 (a) [ The counterex- 
ample (Section |3.3[ ) uses the second column below for the input parameters. 
SFminimizer(6,10) 
xval = 

Columns 1 through 10 



2.1907 


0. 


.0867 


1. 


.0279 2. 


.4529 


0.8294 1.3029 1.0958 0.8502 1.2136 0.3308 


2.1572 


0. 


.4353 


1. 


.9303 0. 


.3841 


-0.2788 0.3006 - 


0.2769 2.8797 3.0782 3.1699 


0.0079 


3. 


.2225 


3. 


.2107 1. 


.8739 


3 


.1482 


3 


.6489 


1, 


.3900 


0.5959 0.1228 1.8107 


0.5496 


2. 


.6633 


0. 


.2000 1. 


.1955 


0, 


.5784 


2, 


.0316 


1, 


.6273 


2.0583 3.1341 1.3276 


1.8511 


2. 


.9336 


1. 


.0424 2. 


.6615 


2, 


.3439 


1, 


.2207 


2, 


.4580 


4.3850 0.1877 -0.1146 


5.2200 


1. 


.2107 


7. 


.0598 3. 


.4826 


0, 


.3696 


0, 


.5633 


5, 


.1432 


4.2546 3.5523 5.5103 


2.7090 


4. 


.1128 


4. 


.1941 2. 


.3445 


2, 


.4291 


2, 


.0589 


0, 


.4295 


2.9467 1.0576 5.9281 


0.1954 


2. 


.0119 


2. 


.2644 7. 


.0575 


2, 


.1168 


2, 


.3027 


1, 


.9744 


-0.2264 2.5206 2.4902 


0.5312 


7. 


.4240 


2. 


.2923 5. 


.1205 


1, 


.2132 


4, 


.0211 


3 


.4493 


1.6729 1.7930 5.7145 


1.3429 


0. 


.8465 


1. 


.4447 4. 


.4661 


4, 


.0120 


1, 


.8860 


0, 


.0947 


3.9222 5.1409 1.2114 


0.5200 


0. 


.2969 


0. 


.7942 0. 


.7334 


0, 


.2149 


0, 


.6042 


1, 


.3269 


0.9360 0.1185 0.0320 


0.0054 


0. 


.0633 


0. 


.2813 0. 


.2543 


0, 


.6229 


0, 


.4588 


0, 


.0165 


0.1667 0.5815 0.5970 



fval = 

Columns 1 through 7 

0.0000 0.0000 1.0000 0.2180 0.0001 1.0000 1.0000 
Columns 8 through 10 
1.0000 0.0001 0.0000 

% The corresponding function values of S and F are given: 
S value = l.Oe - 03* 
-0.3861 -0.0970 0.1462 

Fvalue =2.2329e - 05 
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